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PREFACE 


Shock  test  planning  for  the  Navy  surface  ships  prompted  new  interest  in  the 
phenomena  of  Underwater  Cavitation  as  a  damage  augmenting  process.  Results 
have  been  transmitted  via  letter  reports  during  the  planning  process.  This  is 
the  formal  Final  Report. 


The  following  tasks  were  carried  out: 

1.  Review  past  theory  and  experimental  data; 

2.  Extend  Bulk  Cavitation  theory  for  secondary  pressures  owing 
to  momentum  accretion  and  afterflow  on  bulk  cavitation 
region's  closure  front; 

3.  Extend  theory  to  include  hull  surface's  inertia  and  com¬ 
pliance; 

4.  Assist  in  design  of  bulk  cavitation  experiments;  and 
monitor  DNA's  other  cavitation  efforts. 

Predictions  are  made  for  possible  effects  of  Bulk  Cavitation  in  a  typical 
shock  test. 


Technical  monitors  have  been  Lt.  Robert  Elsbernd,  USN  and  Lt .  Dennis  Sobota, 

USN . 


The  Underwater  Explosives  Research  Division  of  David  Taylor  Research  & 
Development  Center  has  considerable  background  information  on  Bulk  Cavitation 
and  Hull  Plate  Loading.  The  writer  is  grateful  for  recent  help  provided  by  that 
Division's  John  Wise,  Robert  Walker  and  John  Gordon,  and  also  for  early  help  by 
W.  W.  Murray  and  Heinrich  Schauer. 
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SHOCK  INDDCED  CAVITATION 


SECTION  1 
INTRODUCTION 

Daring  shock  test  planning  for  the  Navy  surface  ships,  there  was  renewed  interest  in 
the  phenomenon  of  bulk  cavitation  as  a  damage  augmenting  process.  Results  of  the  subject 
contract  were  transmitted  via  letter  reports  during  the  planning  process.  This  is  the 
formal  Final  Report. 

The  following  tasks  were  carried  out: 

1.  Review  past  theory  and  experimental  data; 

2.  Extend  Bulk  Cavitation  theory  for  secondary  pressures  owing  to  mo¬ 
mentum  accretion  and  afterflow  on  bulk  cavitation  region's  closure 
front; 

3.  Extend  theory  to  include  hull  surface's  inertia  and  compliance; 

4.  Assist  in  design  of  bulk  cavitation  experiments;  and  monitor 
DNA's  other  cavitation  efforts. 

Predictions  are  made  for  possible  effects  of  Bulk  Cavitation  in  a  typical  shock  test. 

Bulk  Cavitation  theory  is  provided  in  Section  2.  The  Closure  process,  with  conse¬ 
quent  water  hammer,  is  described  in  Section  3. 

The  REFfckENCES  cite  the  history  of  Bulk  Cavitation,  Local  Cavitation  and  hull  plate 
loading.  Past  cavitation  experiments  are  discussed  in  Section  4.  Laboratory  experi¬ 
ments,  using  an  exploding  wire  to  simulate  a  nuclear-like  steam  bubble,  are  also  dis¬ 
cussed. 

Work  to  date  by  all  has  involved  the  effects  of  the  direct  shock  wave  only  —  unfet¬ 
tered  by  secondary  shock  waves  owing  to  bottom  reflection.  Bottom  reflections,  see  Sec¬ 
tion  5,  can  play  a  substantial  role.  Further  investigation  is  required. 

Fundamentals  of  hull  plate  loading  associated  with  local  cavitation  are  discussed  in 
Section  6.  Quantitative  detail  involves  diffraction  analysis,  not  yet  achievable  in 
closed  form  analysis.  Simplified,  closed  form  analyses  are  discussed.  Computer  codes 
(including  diffraction  phenomena)  and  results  are  discussed. 

The  Appendix  lists  a  Z100  Basic  code  formulated  to  determine,  for  a  spherical  ex¬ 
plosion,  the  extent  of  Bulk  Cavitation  as  well  as  the  depth  of  final  closure  and  intensity 
of  water  hanmer  shocks  generated. 
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SECTION  2 


BULK  CAVITATION 


The  phenomena  of  Bulk  Cavitation  and  Closure  have  been  described  in  work  dating  back 


.  i  qj  3  1  5  ^  ^  2  7  4  J 

to  1943  *  *  * 


For  a  brief  description  of  the  phenomena  we  follow  reference  15. 


Consider  a  compressive  underwater  shockwave  moving  toward  a  free  water  surface  as 
shown  in  figure  la.  We  consider  a  shock  wave  with  initial  pressure  p0  followed  by  an 
exponential  decay  of  duration  tw-  At  distances  from  the  explosive  source  which  are  of 
practical  interest,  the  velocity  of  the  shockwave  is  essentially  the  speed  of  sound,  c, 
in  water.  Hence,  the  wavelength  X  of  the  incident  shockwave  is 


X  =  ctw 


(1) 


2.1  UPPER  CAVITATION  DEPTH. 

Just  before  shock  front  arrival  at  the  surface,  the  total  absolute  pressure  is  as 
depicted  in  figure  lb.  When  the  compressive  front  arrives  at  the  surface,  boundary  con¬ 
ditions  require  that  a  rarefactive  wave  be  reflected.  Just  after  reflection  from  the 
surface,  the  total  pressure  is  as  depicted  in  figure  2.  We  note  that  as  the  rarefactive 
front,  of  strength  -p0,  moves  down  into  the  water:  (1)  the  pressure  beneath  the  rare¬ 
faction  front  is  that  owing  to  the  waning  portions  of  the  incoming  wave;  (2)  above  the 
rarefaction  front,  the  total  pressure  remains  pa  (atmospheric  pressure)  at  the  surface 
and  progressively  falls  as  the  rarefaction  front  is  approached.  Below  some  point,  shown 
as  Z0,  the  total  absolute  pressure  is  negative  —  the  water  goes  into  tension. 

Since  seawater  cannot  withstand  appreciable  tension,  the  water  ruptures;  the  pres¬ 
sure  in  the  ruptured  water  falls  to  vapor  pressure  —  zero  for  all  practical  purposes. 

All  types  of  incident  shock  waves  do  not  necessarily  result  in  cavitation.  If  the 
peak  pressure  is  sufficiently  large  and  the  decay  wavelength  X  is  exceedingly  long,  hy¬ 
drostatic  pressure  will  rise  faster  than  the  incident  wave  decays,  and  total  pressure 
will  monotonically  increase  with  depth  —  as  indicated  by  figure  3's  curve  III.  A  shor- 
ter  wavelength  may  result  in  pressure  lowering  as  shown  in  curve  II,  but  zero  pressure  is 
never  reached.  Most  practicable  explosions,  including  nuclear,  have  pressure/wavelength 
combinations  (by  the  time  the  free  water  surface  is  encountered)  such  that  zero  pressure 
is  achieved  and  bulk  cavitation  commences  at  some  depth  Z0> 

Figure  4  depicts  an  incident  shockwave  rising  vertically  to  the  free  water  surface 
Zs.  which  is  initially  at  zero  depth.  If  the  incident  shock  has  a  particle  velocity  of 
u0,  upon  reflection  the  surface  Zs  initially  moves  upward  with  velocity  2uo*  When  the 
reflected  wave  reaches  the  onset  depth  Z0,  the  water  ruptures;  total  pressure  at  the 
depth  is  substantially  zero  —  the  vapor  pressure  of  water  at  ambient  temperature. 
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PRESSURE  P  WATER  SURFACE 


Figure  1.  (a)  Compressive  underwater  shockwave  moving  toward  a  free  surface.  Peak 

pressure  pQ;  decay  time  constant  tw.  (b)  Total  absolute  pressure  just 
before  shock  arrival  at  free  surface. 


WANING  PORTION  OF  INCIDENT 
SHOCKWAVE 
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OF  INCOMING 


Figure  4.  Notion  of  water  particle*  at  various  depths,  ranging  item  z,  (the  water  sur¬ 
face)  to  *0  (the  depth  of  cavitation  onset),  zc  (depth  of  ultiaate  cavitation 
closure)  and  zj  (terminal  depth  beyond  which  there  is  no  cavitation).  The 
particle  originally  at  depth  zi  is  typical  of  the  water  particle  which  is 
accreted  to  the  lower  boundary  zj,  of  the  spall.  The  water  particle  of  depth 
*2  follows  a  ballistic  trajectory  till  it  falls  back  to  its  original  depth. 
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2.2  ACCRETION  TO  SPALL 


At  that  instant  there  is  a  liquid  layer  of  water  —  a  spall  of  thickness  (Zs-Z0)  — 
moving  upward  with  almost  the  velocity  2u0.  The  decelerating  forces  are  two:  (1)  boun¬ 
dary  forces  —  atmospheric  pressure  at  the  upper  surface  and  negligible  vapor  pressure  at 
the  lower  surface;  and  (2)  body  forces  —  gravitational  force  acting  throughout  the  li¬ 
quid  layer. 

As  shown  in  figure  4,  the  spall's  lower  surface  is  denoted  by  Zj (Zs-Zl)  —  the 
spall's  thickness  —  is  initially  equal  to  (Zs-Z0) •  We  shall  see  that  this  thickness 
grows.  As  it  does  the  added  deceleration  owing  to  atmospheric  pressure  at  the  upper  sur¬ 
face  becomes  relatively  less  compared  with  the  body  force  of  gravity.  Hence  the  spall 
deceleration  will  asymptotically  approach  lg. 

The  particle  velocity  u  at  any  point  of  the  incoming  shockwave  is,  since  the  acoustic 
approximation  is  adequate 

u  =  p/pc  ,  (2) 


where  p  is  the  water  density,  and  c  is  the  speed  of  sound.  As  the  reflected  rarefaction 
front  moves  into  water  deeper  than  that  at  Z0,  it  encounters  the  waning  portions  of  the 
incoming  wave  —  where  the  particle  velocity  has  fallen  to  u,  with  u  described  by  equation 
2.  Accordingly,  the  rarefaction  front  leaves  behind  it  cavitated  water  whose  initial 
upward  velocity,  or  launch  velocity,  is  2u.  Hence  the  launch  velocity  continuously  falls 
as  the  rarefaction  penetrates  the  waning  portions  of  the  incoming  wave. 

Since  the  cavitated  region  is  entirely  at  vapor  pressure,  the  deceleration  of  the 
launched,  cavitated  water  is  not  affected  by  pressure  forces,  only  gravitational  forces. 
Hence  each  cavitated  water  particle  follows  a  ballistic  trajectory  —  with  exactly  one  g 
of  deceleration. 

We  recall  that  the  liquid  spall  decelerates  asymptotically  toward  lg  —  at  all  times 
it  decelerates  faster  than  lg.  Hence  the  liquid  spall  is  continuously  decelerating  into 
underlying  cavitated  water.  The  spall  thickness  grows  owing  to  accretion  of  the  under¬ 
lying  cavitated  water.  Figure  4  depicts  that  the  spall's  lower  surface  ZL  contin¬ 
uously  colliding  with  underlying  spalled  water.  Spall  growth  continues  in  this  fashion 
until  it  accretes  the  cavitated  particle  which  had  been  launched  from  depth  Zc,  from  what 
we  shall  see  is  the  cavitation  closure  depth. 


2.3  LOWER  CAVITATION  DEPTH. 

The  reflected  rarefaction  front,  as  shown  in  figure  4,  continually  produces  cavi¬ 
tation  until  some  lower  depth  Zj,  the  cavitation  termination  depth.  Beyond  that  depth 
the  hydrostatic  pressure  is  too  large  to  allow  the  total  absolute  pressure  to  reach  zero. 
The  water  below  Zj  remains  liquid  at  all  times. 

The  cavitated  water  just  above  Zj  has  a  very  small  launch  velocity  and  a  very  brief 
ballistic  trajectory.  It  soon  falls  back  on  the  underlying  liquid  water.  Analogous  to 
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behavior  of  the  accreting  spall,  the  under lying  liquid  surface  moves  upward  '*in*  to  cum¬ 
ulative  accretion  of  water  that  had  been  in  ballistic  trajectory  in  the  cavitated  region. 


2 .4  WATER  HAMMER, 

The  lower  cavitation  boundary  rises  while  the  spall's  lower  surface  Z l  falls  —  until 
eventually  there  is  no  cavitated  region  left.  He  particle  velocity  on  the  lower  cavi¬ 
tation  boundary  is,  to  a  first  approximation,  zero.  The  particle  velocity  on  the  upper 
cavitation  boundary  is  effectively  the  velocity  of  the  decelerating  spall.  Hence  at  time 
tc  there  is  a  water  hammer  as  the  spall  impacts  the  quiescent  lower  boundary.  (Actually 
the  particle  velocity  of  the  lower  boundary  is  not  exactly  zero.  The  computer  program 
listed  in  the  Appendix  takes  the  actual  velocity  into  account.) 

Figure  5  is  identical  to  figure  4  —  except  that  the  velocity  of  the  incoming  and 
reflected  fronts  is  assumed  to  be  infinite;  i, e.,  the  duration  of  the  cavitation- 
producing  process  is  assumed  negligible  compared  with  spall  flight  duration  and  the 
cavitation  closure  and  water  hammer  processes.  This  simplification  has  enabled  a  closed 
form  solution  1S,4‘  Qf  the  bulk  cavitation  process.  Results  compare  well  with  experiment 

4  3  4  7 
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As  shown  in  figure  4,  time  t  is  measured  from  the  moment  the  incoming  shockwave  reach¬ 
es  the  water  surface.  Depth  Z  is  measured  from  the  initial  position  of  the  free  surface. 

The  figure  shows  the  Z-t  trajectories  of  the  incoming  shock  front  and  also  the  re¬ 
flected  rarefaction  front.  These  waves  move  at  velocity  c,  the  speed  of  sound  in  water  - 
-  about  1,500  m/sec.  Cavitation  conditions  exist  when  the  reflected  wave  reaches  the 
onset  depth  ZQ.  Below  that  depth  the  water  vaporizes  when  the  rarefaction  front  passes. 

Ultimately  the  upper,  spalled  layer  will  hammer  the  underlying  liquid  water  beneath 
the  boundary  Zg;  the  dynamics  of  the  process  produce  secondary  underwater  pressure  waves 
owing  to  this  progressive  water  haomer. 


2.5  TRAPPED  ENERGY. 

Depending  on  explosive  yield  and  depth  of  burst,  a  large  fraction  of  the  total  inci¬ 
dent  shockwave  energy  can  be  stored  temporarily  as  kinetic  energy  and/or  gravitational 
potential  energy  (of  the  overlying  spall  and  also  of  the  cavitated  particles)  during  the 
bulk  cavitation  process;  and  the  stored  energy  is  re-emitted  later  during  the  water 
hammer  process  when  the  cavitation  closes. 

Bulk  Cavitation's  late-time  secondary  waves  have  produced  extensive  damage  to  ships 
during  tests  at  HARDTACK  Unbrella  and  in  subsequent  HE  tests.  This  prompted  earlier  in¬ 
vestigations.  The  program  for  the  ARKANSAS  shock  tests  prompted  DNA's  renewed  interest 
in  Bulk  Cavitation. 
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Figure  5.  Water  particle  trajectories  under  the  approxiaation  that  the  transit  tiae  of 
the  rupturing  rarefaction  front  is  negligible  canpared  with  the  trajectory 
tiaes  of  the  cavitated  water. 


2.6  SPHERICAL  EXPLOSIONS. 

Figure  6  shows  the  more  practical  situation  of  a  spherical  explosion.  The  incident 
shock  first  arrives  at  the  surface  directly  above  the  explosive,  at  surface-zero  (SZ). 
Launch  velocity  of  the  spall  is  maximum  there:  it  takes  maximum  time  for  the  spall  to 
complete  its  trajectory  before  impacting  with  the  underlying  water. 

As  we  increase  horizontal  distance  from  SZ,  the  incident  shock  arrives  later,  it  is 
weaker  and  its  angle  of  incidence  with  the  surface  results  in  a  decreasing  spall  launch 
velocity.  The  consequence  is  in  figure  7,  which  shows  the  time  of  closure  Tc  as  a  func¬ 
tion  of  horizontal  radius  fi  m  surface-zero.  Also  shown  is  the  time  of  arrival  Tj  of  the 
incident  shockwave.  We  see  that  cavitation  closure  first  occurs  at  radius  at  time 

TC1 . 

Figure  8  shows  the  maximum  surface  excursion  and  the  time  of  maximum  excursion  for 
this  explosion  configuration. 


10 


LAUNCH  VELOCITY 
OF  SPALL 


SHOCK  ARRIVAL  TIME  AT  SURFACE 


EXPLOSIVE  SOURCE 


Figure  6.  Geometry  for  analysis  of  spherical  explosions. 
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Figure  7.  Some  variables  as  a  function  of  horizontal  range  for  a  spherical  explosion. 
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Figure  8.  Sane  variables  u  .  function  of  horizontal  range  for  a  spherical  explosion. 


SECTION  3 


CAVITATION  CLOSURE 


3.1  PLANAR  SHOCKWAVE. 

An  underwater  shockwave,  at  distances  from  the  explosive  source  which  are  of  inter¬ 
est,  behaves  virtually  as  an  acoustic  wave.  In  practice  the  wave  has  a  shock  rise  to  peak 
pressure  pQ,  followed  by  an  exponential- like  decay.  When  such  an  acoustic  wave  is  re¬ 
flected  from  a  rigid  surface  (1)  the  wave’s  momentum  (total  impulse)  is  reversed,  and  (2) 
the  wave's  kinetic  energy  is  conserved.  A  corollary  of  this  is  that  the  wave's  shape 
(magnitude  and  duration)  is  preserved. 

When  an  underwater  shock  wave  is  reflected  from  a  free  surface,  where  negative  pres¬ 
sure  cannot  exist  in  the  medium,  spallation  or  cavitation  occurs.  Die  wave's  energy  be¬ 
comes  stored  in  the  gravitational  energy  of  the  spall,  and  the  energy  is  returned  (with 
considerable  dissipation)  during  the  closure  process.  The  result  is,  eventually,  a  re¬ 
flected  wave  which  is  compressive.  The  wave's  (1)  kinetic  energy  is  not  conserved,  and 
(2)  while  the  magnitude  of  momentum  is  conserved,  the  momentum  is  reversed.  A  corollary 
of  this  is  that  wave  shape  is  not  preserved. 

For  analysis  we  assmne  the  Instant  Cavitation  depicted  in  Section  2's  figure  5  (page 
9);  i.e.,  the  travel  time  of  the  cavitating  wave  front  —  which  moves  at  sound  speed  — 
is  negligible  compared  with  the  duration  of  spallation  and  cavitation.  Ihe  particle  vel¬ 
ocity  u(Z)  at  depth  Z  in  the  cavitated  region  (for  incident  wave  with  shock  rise  to  pres¬ 
sure  p0,  and  exponential  decay  with  wavelength  X)  is  15 

u(Z)  =  -(pa  +  pgZ  +  2p0e~2Z/*)/pc  ,  (3) 


where  pa  is  atmospheric  pressure,  pgZ  is  hydrostatic  pressure  at  depth  Z  and  pc  is  the 
water's  acoustic  impedance. 

In  the  cavitated  region  the  pressure  is  constant  vapor  pressure  of  water.  There  is 
no  pressure  gradient;  hence,  each  cavitated  water  particle  follows  a  ballistic  trajec¬ 
tory.  The  particle's  flight  time  (the  time  until  gravity  reverses  the  particle’s  mo¬ 
mentum  and  it  falls  back  to  its  original  depth  Z)  is 

tf  =  2u(Z)/g  .  (4) 

For  water  at  depths  above  the  closure  point,  the  water  particle  doesn't  get  a  chance  to 
complete  its  ballistic  flight;  the  higher  acceleration  (owing  to  the  added  force  of  at¬ 
mospheric  pressure  on  top)  of  the  overlying  spall  causes  the  particle  to  accrete  to  the 
spall's  underside  before  flight  time  is  completed.  Particles  launched  at  the  closure 
depth  complete  the  ballistic  flight.  Particles  launched  below  the  closure  point  also  do 
not  complete  the  ballistic  flightp  they  accrete  to  the  rising  underlying  water  surface. 
The  front  of  this  accretion-by-fallback  is  the  bulk  cavitation  lower  closure  front. 

The  water  particle  at  closure  depth  Zc  is  the  last  particle  to  complete  its  bal- 
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listic  flight  time.  That  particle  accretes  its  mass  and  momentum  at  time 


tf(Zc)  =  2u(Zc)/g 


(5) 


For  short  wavelength  shock  waves,  the  hydrostatic  pressure  and  the  atmospheric  pressure 
are  negligible  compared  with  the  'pressure'  component  (see  equation  2)  .  Hence  the  flight 
time  is 


tf(Zc)  =  4p0/(pgc)e  2ZC/X. 


(6) 


3.2  SPHERICAL  SHOCKWAVE. 

In  Section  2  it  was  pointed  out  that  there  is  a  radius  (about  surface  zero)  where 
cavitation  closure  and  water  hammer  first  occur.  Figure  9  shows  the  spall  and  the  under¬ 
lying  water  in  the  neighborhood  of  this  radius  X^  0f  first  closure.  Figure  9a  depicts 
the  situation  just  prior  to  first  impact;  9b  depicts  it  immediately  after  first  impact. 
Closure  progresses  inward  with  instantaneous  contact  point  Xj»  and  outward  with  contact 
point  Xj,. 

Immediately  after  first  impact,  X^  and  Xq  travel  approximately  horizontally  with  in¬ 
finite  velocity  and  thereafter  progressively  slow  down. 

The  outward  travelling  point  Xj,  cannot  move  slower  than  the  horizontal  velocity  of 
the  incident  shockwave,  which  travels  at  the  speed  of  sound  in  water.  Thus  the  outward 
travelling  contact  point  monotonically  slows  to  mach  1.  However,  as  it  approaches  mach  1 
the  distance  from  the  explosion  center  and  the  angle  of  incidence  with  the  surface  are 
such  that  the  energy  stored  in  the  cavitation  process  is  negligible,  and  the  consequent 
water  hanner  at  closure  is  negligible. 

The  inward  travelling  closure  point  X^  is  not  limited  by  the  speed  of  sound  in  water. 
It  may  fall  to  and  below  mach  1.  Distance  and  angle  of  incidence  considerations  permit 
appreciable  stored  energy  and  final  water  hammer.  Indeed,  the  horizontal  radius  X^  at 
which  the  contact  point  is  travelling  at  mach  1  is  particularly  hazardous  since  the  ham¬ 
mer  pressure  theoretically  increases  without  bound  at  that  radius15. 

In  simple  analyses  it  is  assumed  that  the  horizontal  motion  of  cavitated  water  par- 
tides  is  negligible  compared  with  the  vertical  motion.  Figure  10  depicts  the  mach  wave 
generated  by  the  supersonically  (in  water)  travelling  closure  point.  Figure  10  shows  the 
vertical  velocities  of  the  spall  (us)  and  the  underlying  water  (u^) .  Figure  10's  equa¬ 
tion  4  shows  the  particle  velocity  normal  to  the  mach  front  and  to  the  reflected  mach 
front.  Equation  6  shows  that  the  hasner  pressure  Pjj  increases  without  bound  as  the  clos¬ 
ure  mach  number  M  slows  down  and  approaches  unity.  Even  so,  equations  9  and  10  show  that 
the  hanmer  impulse  IQ  below  closure  depth  Zc  is  bounded.  Above  closure  depth  the  hammer 
impulse  I  is  proportional  to  the  depth  Z  below  the  water  surface. 

If  a  vessel's  structure  is  pressure  sensitive,  it  is  at  risk  of  self  damage  if  it 
launches  an  underwater  depth  charge  so  as  to  place  itself  at  the  mach  1  distance.  A 
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Fi for*  9.  Spall  and  underlying  water  In  the  neighborhood  of  radios  of  first 

closure,  X^,  where  water  banner  first  occurs. 
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Figure  10.  Mach  wave  generated  by  supersonically  travelling  locus  of  closure. 


target  vessel  is  most  vulnerable  if  placed  at  that  location.  There  is  less  hazard  if  the 
vessel's  structure  is  impulse  sensitive. 

Figure  11  graphs  the  closure  mach  nun  be  r  and  the  time  of  closure  for  a  40,000  LB.  H1X 
charge.  The  right  ordinate  is  negative  for  the  inward  travelling  mach  front,  and 
positive  for  the  outward  front. 

Figure  12  graphs  the  water  haunter's  duration,  pressure  and  impulse  for  the  same  ex¬ 
plosion  configuration. 

Figure  13  shows  the  extent  of  bulk  cavitation  and  also  the  locus  of  closure.  The 
data  for  figure  13  resulted  from  the  computer  code  given  in  the  Appendix.  The  Appendix* 
code  is  an  embellishment  of  the  fortran  code  developed  in  reference  4.  As  shown  in  the 
Appendix,  the  code  was  modified  in  reference  9  to  account  for  bottom  layer  acceleration, 
and  modified  for  the  present  report  to  account  for  afterflow  (10). 
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Figure  11.  Closure  aach  msnber  and  time  of  closure  for  spherical  explosion. 
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Figure  13.  Extent  of  bulk  cavitation  and  locus  of  closure  for  spherical  explosion. 


SECTION  4 


EXPERIMENTAL  DATA 


4.1  FIELD  TESTS. 

There  have  been  three  field  tests  specifically  designed  to  investigate  bulk  cavi¬ 
tation.  The  first  was  conducted  at  UERD  in  1962  and  reported  by  Walker  and  Gorden  (ref 
47).  Theoretical  interpretation  of  its  data  has  been  good  —  until  arrival  of  the  first 
bottom  reflections.  Thereafter  data  seems  to  be  the  superposition  of  confusing  rever¬ 
berations. 

The  second  set  consisted  of  the  NOL  Tests  at  Mono  Lake  conducted  in  1969  (ref  43). 
Again,  they  produced  data  agreeing  reasonably  with  theory  during  early  times.  The  Mono 
Lake  tests  placed  floats  with  a  seismic-suspension  integrating  accelerometer  to  measure 
surface  velocity.  Figure  14  shows  typical  results.  The  integrating  accelerometer  dis¬ 
played  excessive  integration  drift  about  70  msec  after  onset  of  motion. 

Data  from  the  Mono  Lake  tests  must  be  used  judiciously  owing  to  unusual . sound  prop¬ 
agation  behavior  associated  with  1)  high  salt  concentration;  2)  substantial  concentra¬ 
tion  of  brine  shrimp. 

The  third  and  more  recent  tests  were  conducted  by  NRL  and  UHID.  They  included  a 
more  recently  developed  capability  to  place  PV  gages  under  water  —  providing  useful  cor¬ 
relations  between  pressure  and  particle  velocity  throughout  the  field. 

Bulk  Cavitation  theories  do  not  accomnodate  bottom  reflection  (discussed  in  Section 
5).  If  there  is  to  be  a  full  understanding  of  Bulk  Cavitation  it  will  be  necessary  to 
have  data  from  at  least  one  test  in  deep  water,  where  measurements  are  not  confounded 
early  on  by  bottom  reflections.  Enough  depth  is  needed  to  make  sure  that  a  complete 
cycle  of  bulk  cavitation  —  incident  pressure,  generation  of  cavitation,  closure,  and 
measurement  of  secondary  waves  —  is  completed  before  arrival  of  the  bottom  reflection. 


4.2  LABORATORY  TESTS. 

Laboratory  tests  were  carried  out  in  1969  using  an  exploding  wire  as  the  explosive 
source  3*.  HE  tests  produce  a  gas  bubble  consisting  of  HE  combustion  products,  whereas 
the  exploding  wire  is  nuclear-like  in  that  it  generates  a  steam  bubble. 

The  facility's  original  use  was  to  gather  information  on  bulk  cavitation.  It  was 
later  put  to  use  to  study  (1)  shockwave/ hull-plate  interaction  with  and  without  local 
cavitation,  and  (2)  surface  wave  generation. 

The  facility  is  shown  in  figure  15.  The  test  vessel  is  a  steel  cylinder  6  feet  in 
diameter,  with  near  hemispherical  ends  to  permit  vacuum  or  pressurization  to  30  psig. 
The  interior  is  lined  with  anechoic  material. 

The  principal  observation  windows  are  diametrically  opposite,  12  inches  in  diameter 
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Figure  14.  Spell  surface  velocity  in  bulk  cavitation  field  test  at  Kmo  Lake.  1969. 
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and  consist  of  schlierin-free  glass.  One  such  window  is  shown  in  figure  16. 

Depending  on  application,  either  schlierin  or  shadowgraph  photos  are  taken.  Light 
sources  available  are  (a)  a  1  psec  spark,  (b)  a  1  msec  xenon  flashlamp,  and  (c)  a  steady 
high  pressure  mercury  arc  lamp. 

The  cameras  for  schlierin  work  are  a  single  frame  4  X  5  camera  as  well  as  a  drum 
camera  taking  224  images  with  frame  separation  of  40  nsec  and  shutter  speed  of  1  psec. 
Figure  17  shows  a  typical  sequence  of  drum  camera  photos. 

For  longer  sequence  durations  a  35  m  Fastax  camera  was  used  —  with  frame  rates  up  to 
3000/sec.  A  16  mm  Bolex  was  also  available  for  very  long  events. 

Miniature  piezoelectric  hydrophones  are  used  to  gather  simultaneous  pressure-time 
information. 

Energy  for  the  explosion  was  provided  by  a  20-KV,  15-pF,  5-nH  capacitor.  The  ex¬ 
ploding  wire  can  be  positioned  anywhere  in  the  tank  —  above  or  below  the  water  surface. 
A  2  inn  length  of  40AWG  nichrome  was  typically  used  —  virtually  a  point  source  for  most 
appl ications. 

The  'nuclear'  yield  of  the  exploding  wire  was  of  the  order  of  0.04  gm  of  TNT.  Data 
analysis  showed  scaling  laws  applied  favorably  over  a  very  wide  range  —  for  full  scale 
nuclear  (CROSSROADS-Baker,  at  the  same  scaled  depth  of  burst  in  shallow  water),  large  HE 
(to  104  lbs  of  tnT)  and  small  HE  charges  (to  10-4  ibs  of  TNT) .** 

With  the  exploding  wire  as  a  simulator  for  nuclear  explosions,  the  schlierin  optics 
show  the  progress  of  direct  and  reflected  shock  waves  in  the  water  and  the  air  and  permit 
observation  of  bottom  effects  and  surface  effects.  Data  connecting  the  generation  of 
surface  waves,  plumes  and  air  blast  have  been  obtained  at  low  cost.  A  variety  of  re¬ 
flecting  bottoms  —  in  material  and  contour  —  can  be  studied. 
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SECTION  5 


BOTTOM  REFLECTION 


5.1  THEORY. 

Officer35  js  a  good  source  of  information  on  acoustic  bottom  reflection.  Figures  18 
through  21  are  interesting  excerpts.  Figure  18  shows  the  geometry  for  explosive  source, 
bottom  and  receiver,  p  is  the  density  of  the  medium  and  c  is  the  speed  of  sound.  For 
incidence  angles  less  than  the  critical,  8C,  the  reflection  coefficient  is  less  than  100 
percent,  i.e.,  there  is  some  energy  transmitted  into  the  bottom.  Beyond  the  critical 
angle  there  is  total  reflection. 

Fig.  19  shows  an  example  of  the  reflectivity  coefficient  and  phase  change  a  as  a 
function  of  incidence  angle  0  —  for  sinusoidal  sound  trains.  The  curves  are  typical 
for  the  usual  bottom  condition  where  C2  >  cj  and  P2C2  >  pjcj  —  where  cj  and  C2  are  the 
sound  speeds  in  water  and  in  the  bottom  respectively,  and  p^ci  an<j  p2c2  are  the  respective 
acoustic  impedances  in  these  media.  In  such  cases  reflectivity  reaches  100  percent  at 
some  critical  incidence  angle  0C;  for  greater  angles  the  reflectivity  remains  at  100 
percent  while  the  phase  a  of  the  wave  train  begins  to  shift  —  reaching  180°  phase  shift 
at  grazing  incidence,  8  =  90°. 

Such  phase  shifting  causes  distortion  of  an  incident  wave.  For  the  practical  case 
of  a  wave  with  a  shock  rise  followed  by  an  exponential-like  decay,  reference  2  describes 
the  distortion  for  angles  of  incidence  beyond  the  critical,  as  shown  in  figure  20.  (N.B, 

the  angle  a  shown  in  figure  20  is  the  phase  change  angle,  which  is  related  to  the  angle  of 
incidence  8  as  typically  shown  in  figure  19.)  The  phase  angle  remains  zero  —  and  the 
wave  is  undistorted  —  from  normal  incidence  all  the  way  to  the  critical  incidence  angle. 
Thereafter  the  phase  begins  to  shift  and  the  wave  distortion  is  as  shown  in  figure  20. 

When  the  phase  shift  is  larger  than  about  90°  for  this  explosion-like  wave  shape,  one 
observes  a  compressive  precursor  (for  times  t  <  0) ,  and  a  rarefaction  follow-up  —  and 
these  are  for  for  bottom  incidence  angles  generally  encountered  in  underwater  explosion 
tests.  This  may  explain  the  'strange'  bottom  reflection  rarefactions  noted  in  Walker  and 
Gordon's  1962  tests. 

Fr iedlander1*  provides  theoretical  treatment  of  exp 1 os ion- like  pulses.  Figure  21 
shows  his  treatment  of  a  spherical  pulse  reflecting  from  a  bottom.  Reflection  front  B  is 
for  less  than  critical  angle;  front  L  is  at  the  critical  angle.  Front  LEN  is  beyond  the 
critical  angle.  For  example,  along  the  line  ED,  we  have  (1)  a  gradual  rise  compressive 
precursor  (compare  Figure  20)  between  E  and  the  'reflected  front'  and  (2)  between  the 
reflected  front  and  D  there  is  a  main  wave,  compressive  or  rarefactive  (compare  Fig.  20). 
LEN  is  an  unequivocal  precursor  for  that  portion  beyond  the  incident,  near  N. 
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Figure  18.  Bottan  reflection  geometry,  showing  explosive  source,  bottom  and  receiver. 


Reflect 


A  //5> 

a4-  /&■? 

\/Qr 

'  1 

i 

0.2- 


15*  30*  45* 

Angle  of  Inc 


Figure  19.  Phase  change  a  and  reflectivity  coef 
9  —  for  sinusoidal  sound  trains. 


Figure  21.  Friedlander' s  treatment  of  spherical  pulse  reflecting  from  bottom. 


5.2  SHOCK  TESTING  CONSIDERATIONS 


Bottom  reflections  can  play  a  significant  role  in  shock  testing.  Lacking  adequate 
knowledge  of  bottom  conditions,  one  must  meantime  consider  what  can  happen  with  100  per¬ 
cent  reflection,  coupled  with  bulk  cavitation  effects. 

Among  the  test  configurations  considered  was  the  following: 

Explosive  weight:  40,000  lbs  HBX 

Explosive  depth:  200  feet 

Water  depth:  700  feet 

Horizontal  Range:  722  feet 

Using  the  analysis  of  reference  15  (embodied  in  the  program  listed  in  the  Appendix), 
at  a  horizontal  range  of  722  feet  the  bottom  reflected  wave  would  appear  to  cause  a  second 
round  of  bulk  cavitation  (beginning  just  after  closure  of  the  first  cavitation).  The 
impact  velocity  —  and  consequent  closure  hammer  pressure  —  is  larger  than  that  of  the 
first  cavitation  closure.  Further,  the  geometry  is  such  that  tV  second  cavitation  clos¬ 
ure  first  occurs  at  a  horizontal  range  of  1800  feet;  however,  the  mach  1  closure  point  — 
where  peak  pressure  increases  without  bound  —  occurs  at  a  horizontal  range  of  720  feet. 

At  a  horizontal  range  of  722  feet,  the  surface  incidence  angle  of  the  main  shock  is 
about  75°.  The  surface  incidence  angle  of  the  bottom  reflection  is  about  30°.  Even 
though  the  reflected  pressure  is  only  half  (assuming  unity  bottom  reflectivity)  of  the 
main  pressure,  the  change  in  incident  angle  more  than  compensates,  generating  higher 
spall  velocity  and  consequently  higher  closure  hammer  velocity. 

Figure  22  shows  the  expected  pressure  history  at  a  depth  of  21  feet  and  horizontal 
range  of  722  feet.  The  second  cavitation  closure  pressure  is  theoretically  infinite; 
however,  the  impulse  is  finite,  estimated  to  be  about  3  psi-seconds  (compared  with  a 
value  of  0.92  for  the  impulse  owing  to  closure  of  the  first  cavitation). 

Without  adequate  knowledge  of  bottom  conditions  in  proposed  shock  testing  it  is  not 
possible  to  predict  whether  or  not  there  might  be  a  suiprisingly  large  impulse  delivered 
owing  to  interaction  of  bottom  reflection  and  bulk  cavitation. 


CHARGE  40,000  LBS,  HBX 
BURST  DEPTH  200  FEET 
BOTTOM  DEPTH  700  FEET 
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Figure  22.  Possible  expected  pressure  history  for  shock  test  of  ship. 


SECTION  6 


HULL  PLATE  LOADING 

A  prime  concern  for  underwater  shockwave /ship  interaction  is  the  effects  on  equip¬ 
ment  within  a  ship.  Secondly,  observed  hull  plate  buckling  is  enough  to  cause  flow  noise 
that  compromises  a  ship's  sonar  ability.  The  mechanisms  of  hull  plate  buckling  are  not 
adequately  understood. 

Analytical  solutions  to  the  problem  were  first  sought  during  WWII.*5,*7,31,45  Tests 
and  analysis  pointed  to  local  cavitation  playing  a  major  role.  The  underwater  Explosives 
Research  Division  of  NSRDC  conducted  a  number  of  HE  tests  during  the  1960s  to  study  the 
mechanisms  of  hull  damage  from  underwater  explosions.  At  that  time  it  was  realized  that 
there  were  three  intertwined  phenomena  contributing  to  hull  plate  damage:  (1)  initial 
shockwave/hull-plate  interaction,  (2)  local  (and  perhaps  also  bulk)  cavitation  and  clos¬ 
ure,  and  (3)  diffraction. 

The  complexity  requires  machine  numerical  methods.  Even  so,  a  semi-quantitative  a- 
nalysis  —  not  involving  diffraction  —  enabled  interpretation  of  the  UffiD  data.1* 

Though  no  analytical  solution  is  available,  the  differential  equations  representing 
the  problem  can  be  set  up  and  solved  digitally.  Useful  results  have  been  obtained  by  a 
number  of  workers  for  the  situation  where  the  water  is  always  an  acoustic  medium,  i.e,  at 
depth  where  ambient  pressure  allows  large  negative  acoustic  pressures.4 »* ,17 ,*1,- 
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More  recent  numerical  treatment  includes  provision  for  a  bilinear  acoustic  medium 
interacting  with  simplified  structural  configurations. 5 »** »** »34 »* 7 


Original  insights  into  hull  plate  loading  stylized  the  interaction  to  be  that  of 
shockwave  and  plate.4*  Late  time  effects  in  these  analyses  pointed  toward  local  cavi¬ 
tation  as  a  major  affect2*.  Afterflow  was  also  considered  as  a  late  time  loading  mech- 
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6.1  OBSERVATIONS 

Streak  camera  tests  have  been  conducted  on  a  circular  diaphragm  —  held  on  a  test  hull 
bottom  by  a  heavy  rim  —  and  subjected  to  an  HE  shockwave.24  The  results  are  summarily 
sketched  in  figure  23.  The  first  sketch  shows  the  diaphragm  inmediately  after  it  has 
been  struck  by  the  shock  front  —  it  has  just  begun  to  move.  The  space  between  the  dia¬ 
phragm  and  the  reflected  shock  front  has  been  expanded  in  the  adjacent  box  which  shows  the 
pressure  distribution.  There,  we  see  that  the  reflected  shock  front  is  moving  into  the 
waning  portions  of  the  initial,  incoming  shockwave.  The  sum  of  pressures  in  the  incident 
wave  and  the  reflected  wave  cause  substantially  a  pressure  doubling. 

As  the  diaphragm  accelerates  to  higher  velocity,  the  reflected  wave,  initially  a  com¬ 
pression,  becomes  increasingly  rarefactive.  In  the  pressure  distribution  box,  we  there¬ 
fore  observe  the  rapidly  decreasing  total  pressure.  If  water  could  withstand  tension, 
the  total  pressure  would  follow  the  dotted  line.  However,  it  cannot;  hence,  near  the 
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Figure  23.  Smeary  of  streak  camera  tests  showing  response  of  circular  dia¬ 
phragm  to  incident  shockwave. 
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diaphragm  the  water  cavitates  —  by  forming  vapor  bubbles,  or  by  forming  water  droplets 
in  a  matrix  of  water  vapor.  Whichever  way,  the  pressure  in  the  cavitated  region  is  that 
of  water  vapor,  effectively  zero  for  our  purposes. 

Sketch  2  shows  the  diaphragm  at  a  later  time.  The  flat  central  portion  of  the  dia¬ 
phragm  does  not  yet  know  that  the  edges  have  been  clamped  by  the  heavy  rim;  a  bending 
wave,  travelling  inward,  carries  this  information.  The  streak  camera  record  displays  a 
jerky  motion  as  shown  in  the  second  and  later  sketches.  Theoretical  work  has  not  shown  a 
cause  for  this  jerkiness.  It  is  conjectured  that  a  water  layer,  still  attached  to  the 
diaphragm,  is  vibrating. 

Later  sketches  show  further  motion  and  progress  of  the  bending  wave.  Next  to  last 
sketch  shows  the  bending  wave  reaching  the  center.  The  camera  record  indicates  that  all 
motion  has  momentarily  stopped. 

Then,  at  a  substantially  later  time,  the  diaphragm  is  bulged  outward.  It  appears 
that  this  late  impulse  is  a  water  hammer  that  takes  place  at  cavitation  closure. 

6.2  PRELIMINARY  ANALYSIS. 

Kennard  has  displayed  the  physics  of  shockwave/hull  interaction  involved  quite  well 
by  a  variety  of  analytical  treatments  with  restricting  assumptions.2 7 

Four  characteristic  times  are  used  (see  figure  24b): 


1  .  ty, 

2 .  tm 

3.  td 

4 .  ts 


wave  time  —  the  incident  wave  is  assumed  for  analytical 
purposes  to  have  a  shock  rise  to  pressure  po  followed  by 
exponential  decay  with  time  constant  tw. 

maximum  time  —  required  for  the  hull  structure  being  ana¬ 
lyzed  to  reach  maximum  velocity. 

diffraction  time  —  the  time  required  for  an  underwater 
acoustic  wave  to  propagate  fran  the  structure's  center  to 
its  edge. 

swing  time  —  the  time  required  for  the  structure  to  reach 
maximum  deflection  and  come  to  rest. 


The  simplest  analysis  showing  hull/cavitation  phenomena  is  that  of  an  exponential¬ 
like  shock  wave  striking  a  plate  of  infinite  lateral  extent,  as  shown  in  figure  24a. 2 7 >4 5 


If  the  areal  density  of  the  plate  is  m,  a  damping  time  constant  tp  comes  into  play, 
defined  by 


tp  =  m/pc  .  (7) 

tp  may  be  visualized  as  the  time  required  for  a  sound  wave  to  traverse  a  thickness  of 
water  having  the  same  mass  as  the  plate.27  For  example,  a  one- inch  steel  plate  has  a 
damping  time  of  0.129  msec;  a  rigid  hull  drawing  30  feet  of  water  has  an  effective  damp¬ 
ing  time  of  6  msec  (for  a  wave  shocking  it  from  the  bottom) . 

In  figure  24a  we  see  the  incident  shockwave.  At  any  instant  at  the  plate  the  inci¬ 
dent  plus  reflected  shock  wave  sun  to  a  pressure  doubling,  i.e,  a  pressure  of  2pj. 
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(a)  (b) 


Pi gore  24.  (*)  Exponent ial- like  shockwave  incident  on  plate  of  infinite  lateral  extent, 

(b)  Characteristic  tines  employed  in  analysis  of  shocked  plate. 
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Superposed  on  this  is  a  rarefaction  wave  owing  to  motion  of  the  plate.  As  indicated  in 
figure  24a,  the  magnitude  of  the  rarefaction  is  pc*dz/dt  —  where  p  is  the  water  density 
(65  lbs/ft3  or  2  slugs/ft3),  c  is  the  speed  of  sound  in  water  (5000  fps),  and  z  is  the 
displacement  of  the  plate,  pc  is  the  acoustic  impedance  of  water. 


The  equation  of  motion  for  the  problem  shown  in  figure  24a  is 


m’d^z/dt^  =  2pj  -  pc*dz/dt 


(8) 


Hie  incident  shock  pressure  is  described  by 


Pi  *  P0e-t/tw  • 


(9) 


Wave  time  t^  and  peak  pressure  p0  follow  empirical  scaling  laws3 


tw  =  MW1/3(R/W1/3)®  , 

(10) 

p0  =  K(W1/3/R)k  . 

(11) 

where  R  is  the  radius  from  the  explosive  and  the  empirical  constants  are  shown  in  Table  1 
for  various  explosives. 10 «11»41 

Table  1 


Explosive 

K,  psi 

k 

M,  msec 

m 

HBX-1 

23.980 

1.13 

0.0654 

0.18 

Lithanol 

13,180 

1.13 

0.12 

0.18 

Pentolite 

22,500 

1.13 

0.0969 

0.08 

TNT 

21,600 

1.13 

0.0676 

0.24 

The  ensuing  graphs  of  pressure  and  motion  of  the  shockwave /plate  interaction  are 
plotted  in  terms  of  the  following  dimensionless  variables: 


r  -  tp/tw  , 

(12a) 

^*i  =  Pi/Po  » 

(12b) 

pt  =  Pt/Po  * 

(12c) 

Pr  =  Pt  '  pr  • 

(12d) 

U  =  l/(pQ/pc)  , 

(plate  velocity) 

(12e) 

% 

^  * 
ii 

H 

(12f) 

KE  =  1  PfU.dT  , 

Jo 

(kinetic  energy  of  plate) 

(12g) 

Z  =  z/X  , 

(12h) 

where  shock  wavelength  X  is  related  to  t*  by 
X  =  c.t*  , 


( 12  i) 
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and  c  is  the  speed  of  sound  in  water. 


In  terms  of  these  the  dimensionless  input  shockwave  is  expressed  by 


Pi  =  e_T  , 

(13  a) 

various  responses  are: 

Pt  =  2[e-T/r  -  re~T]/ ( l-r)  . 

(13b) 

Pr  -  Pt  -  Pr  • 

(13c) 

U  =  2[e-T  -  e-T/r]/(i-r)  . 

(13d) 

These  are  plotted  in  figure  25  for  respective  r  values  of  0,  1/4,  1  and  4. 


6.3  LOADING  AT  DEPTH. 

If  plate  loading  takes  place  at  sufficient  depth  —  ambient  pressure  (atmospheric 
plus  hydrostatic)  is  large  —  then  the  acoustic  wave  negative  pressures  are  realizeable: 
the  entire  plate  loading  process  can  be  described  in  terms  of  the  foregoing  graphs. 

The  area  under  the  Pt  curve  is  the  impulse  imparted  to  the  surface  plate.  The  im¬ 
pulse  is  maximum  at  some  time  Tto,  depending  on  the  value  of  the  parameter  r.  Thereafter 
Pt  goes  negative  and  thereby  begins  decelerating  the  surface  plate;  it  begins  to 
withdraw  impulse  from  the  plate.  Indeed  if  we  wait  several  time  constants  —  let  T 
approach  infinity  —  the  acoustic  pressures  will  withdraw  all  impulse  loading  from  the 
plate.  The  plate  will  came  to  rest. 

Figure  26  shows  values  of  Ttc  for  values  of  r. 

Pt>  the  reflected  pressure,  overall  has  a  total  (dimensionless)  impulse  of  zero. 
Hence,  the  total  impulse  of  Pr  is  opposite  to  that  of  Pj.  At  depth,  then,  the  reflected 
pressure  can  radiate  downward  to  infinity,  carrying  away  energy  and  momentum,  ultimately 
leaving  the  plate  motionless. 

Hence  in  a  totally  acoustic  plate  loading,  the  maximum  impulse  (up  to  a  maximum  of  -2 
for  large  r)  is  imparted  to  the  plate  at  time  Tt0.  Thereafter  the  acoustic  waves  proceed 
to  unload  the  plate,  ultimately:  (1)  leaving  it  with  zero  imparted  impulse;  (2) 
delivering  all  the  incident  wave's  impulse  to  the  reflected  wave's  impulse,  which 
radiates  away  to  infinity. 


6.4  LOADING  NEAR  THE  SURFACE. 

If  we  are  not  at  significant  depth  —  i.e.,  a  vessel  at  or  near  the  surface  —  the 
above  described  acoustic  waves  can  be  followed  only  up  until  approximately  Tto>  wben  the 
total  pressure  passes  through  zero.  Thereafter  the  total  acoustic  pressure  Pt  cannot 
continue  into  its  negative  phase;  the  water  rn;;ures  and  cavitation  takes  place. 
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We  can  use  these  acoustic  curves  to  find  the  condition  of  the  plate  at  time  Tt0.  when 
the  plate  has  received  maximum  impulse.  After  cavitation  the  actual  pressure  Pt  remains 
at  zero  (actually  the  vapor  pressure  of  water  —  practically  zero);  no  impulse  is 
subsequently  withdrawn  from  the  plate. 

By  time  T0  the  reflected  wave  Pr  has  become  negative.  This  reflected  wave  is 
travelling  down  into  the  water;  hence,  its  argument  as  a  travelling  wave  is  T-Z 
(dimensionless  time  minus  dimensionless  displacement).  Since  this  reflected  pressure 
wave  Pr(T  -Z)  is  propagating  back  down  into  the  water,  it  moves  into  the  waning  portions 
of  the  incoming  Pj(T  +  Z).  Thus  at  some  place  and  time,  Pr(T-Z)  plus  Pj(T+Z)  becomes 
negative  —  physically  not  allowed  at  shallow  depths  (where  ambient  pressure  —  sum  of 
atmospheric  and  hydrostatic  pressures  —  is  effectively  zero);  the  water  cavitates.  As 
time  progresses  the  boundaries  of  this  growing  cavitated  region  are  detemined  by  the  two 
values  of  Z  that  approximately  satisfy,  at  each  time  T11*15'44 

Pjtt+Z)  +  Pr(T-Z)  =  0  .  (14) 

The  only  portion  of  Pr(T-Z)  that  ultimately  escapes  the  cavitation  process  is  the 
compressive  front,  which  is  seen  in  figure  25  to  end  at  some  time  Tro»  depending  on  the 
value  of  the  parameter  r. 

Beyond  ^ro  figure  25  shows  that  the  additional  impulse  in  the  reflected  wave  Pr(T~Z) 
is  negative.  All  of  this  negative  impulse  —  initially  the  tail  of  Pr(T-Z)  radiating 
down  away  from  the  plate  —  cannot  escape  and  ultimately  is  deposited  in  a  spall  of 
cavitated  water  beneath  the  plate. 

There  are  four  impulses  to  consider  when  the  loading  and  cavitation  processes  have 
been  completed  (but  prior  to  the  beginning  of  the  cavitation  closure  process,  i.e.,  the 
cavitation  closure  phase  of  further  loading  of  the  surface  plate): 

1.  Ij  -  the  total  impulse  of  the  incident  compression  wave  (equal  to 

unity  in  dimensionless  terms); 

2.  If n  -  the  total  impulse  imparted  to  the  surface  plate  by  the  re¬ 

flected  wave; 

3.  Iro  -  the  total  impulse  (which  escapes  to  infinity)  in  the 

compressive  front  of  the  reflected  wave;  and 

4.  Is  -  the  total  impulse  or  momentum  in  the  cavitated  slug  of  water 

immediately  beneath  the  surface  plate. 

The  mathematical  definitions  of  these  are  shown  in  figure  27.  Momentum  conservation 
requires 


li  +  I(o  +  ItQ  +  1$  0  .  (15) 

These  impulses,  as  a  function  of  the  time  constant  ratio  r,  are  graphed  in  Figure  27. 
Also  shown,  in  Figure  26,  are  the  pertinent  times,  Tt0  and  Tro»  as  a  function  of  r. 

r  is  small  if  (1)  the  surface  load  m  is  sufficiently  small  (including  zero,  the 
situation  for  Bulk  Cavitation),  regardless  of  tw  (iie>>  for  mAn 

or  large  explosive 

charges);  or,  (2)  ^  is  large  (corresponding  to  large  scale  HE  and  nuclear). 

For  large  r  cavitation  plays  no  important  role  -  negligible  incident  impulse  is 
deposited  in  a  cavitated  spall  —  cavitation  is  effectively  nonexistent.  Behavior  is 
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Figure  27.  Lnpulses  as  a  function  of  time  constant  ratio  r  (=  tp/tw)  —  for  shocked  plate. 


similar  to  the  well  known  results  for  acoustic  reflection  from  a  rigid  wall:  (1)  the 
reflected  wave  has  an  impulse  equal  in  magnitude,  opposite  in  direction,  to  the  incident 
impulse;  and  (2)  the  wall  receives  an  impulse  equal  to  twice  the  incident  impulse. 

6 .4 .1  QUALITATIVE  DESCRIPTION  OF  PLATE  MOTION. 

Figure  28a  shows  a  shock  wave  impinging  perpendicular  to  hull  plating  supported  by 
stringers.  As  an  approximation  we  consider  the  configuration  of  figure  28b.  Instead  of 
a  diaphragm  with  bending  wave,  we  connect  the  light  weight  hull  plating  by  means  of  a 
spring  to  the  heavy,  rigid  hull  frame.  We  look  at  the  beginning  motion,  while  the  effect 
of  the  spring  is  small. 

Figure  28b  implies  that  the  fraction  f  of  the  water/hull  interface  surface  has 
surface  load  m;  fraction  (1-f)  has  surface  load  M.  The  simplif ication  we  employ  is  that 
the  total  acoustic  pressure  pt  at  each  time  the  same  in  (1)  the  area  fraction  f,  and 
(2)  the  area  fraction  (1-f).  In  other  words  we  assume  that  it  takes  negligible  time  for 
pressures  to  diffract  between  the  area  fraction  f  and  the  area  fraction  (1-f) . 

An  example  that  was  considered  consists  of  a  40,000  pound  charge  of  HBX-1  at  a 
distance  of  600  to  800  feet  from  a  ship  with  1/2  inch  plating  and  drawing  20  feet  of 
water.  Then  we  have 


Po 

500 

psi 

tw 

2.9 

msec 

pD 

1,280 

lbs/ ft^ 

f 

0.8 

m 

20.8 

lbs/ft2  . 

The  overall  ship  mass  is  pD  where  p  is  the  density  of  seawater  and  D  is  the  ship 
draft.  The  loadings  of  the  partial  areas  f  and  (1-f)  must  collectively  support  this 
mass : 


f.m  +  (1-f) -M  =  pD 


(16) 


so  that 


M  =  (pD  -  f.m) / (1-f)  .  (17) 

Since  the  hull  plate  mass  m  is  generally  very  small  compared  with  the  ship  mass,  we  have 
to  a  good  approximation 

M  =  pD/ (1-f)  .  (18) 

If  the  unsupported  portion  of  the  hull  plating  effectively  makes  up  80  percent  of  the  hull 
area,  then  the  massive  loading  M  is  5  times  the  ship  mass  per  unit  area  pD: 

M  6,400  lbs/ft2  . 
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impinging  on  supported  hull  piste,  (b)  Analytical  approximation  of  hull  plate. 


The  damping  time  constants  for  m  and  M  are 


tpm  2.08  msec 

tpM  640  msec  , 

and  the  dimensionless  ratios  r  are 


rM  220  • 

With  these  values  we  see  from  the  graphs  of  figures  25  through  27  that  the  stiff  hull 
portions  are  initially  absorbing  all  the  initial  impulse;  the  reflected  wave  is  high 
pressure.  Ihe  compliant  plate  portion  rapidly  begins  to  radiate  a  low  pressure  wave  and 
traps  much  of  its  energy  in  the  kinetic  energy  of  spalled  or  cavitated  water. 

As  seen  in  figure  29,  the  high  pressure  reflected  wave  pumps  its  energy  into  the  low 
pressure  region,  exacerbating  the  trapped  kinetic  energy  there.  At  some  late  time  this 
kinetic  energy  must  be  absorbed  by  the  hull  plate. 

Descriptions  of  major  effects  have  been  provided  here  to  the  extent  feasible  without 
explicitly  including  diffraction  time  considerations  in  the  analysis.  Rigorous  answers 
require  machine  numerical  solution. 
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Figure  29.  Beginning  motion  of  supported  hull  plate,  showing  pooping  of  high  pressure  into  cavitated  region. 
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APPENDIX 


BULK  CAVITATION 
ZIOO  Basic  Program 


*******  FORTRAN  VERSION  OF  1 CAVIT '  **** 
ORIGINAL  PROGRAM  BY  COSTANZO  AND  GORDON 
UERD/DINSRDC  (DNA) .  » 

9  OCT,  1980  REVISION  TO  INCLUDE  » 
EFFECTS  OF  BOTTOM  LAYER  ACCELERATION* 
BY  J.  CLARK  (CUA).  * 

23  DEC  1987  REVISION  TO  INCLUDE  * 
UPPER  BOUNDARY  DYNAMICS,  AFTERFLOW  « 
AND  CLOSURE  MACH  NUMBER  * 

BY  V.  CUSHING  (CAIJ.  * 


10  '  CAVIT. 7CN  BASIC  VERSION 

12  GOTO  40 

20  CLS :  LOCATE  4,1 

22  PRINT  "  »««**»*»»»  FORTRAN  VERSION  OF  »CAVIT'  *»***»«" 

24  PRINT  "  ***ORIGINAL  PROGRAM  BY  COSTANZO  AND  GORDON***" 

26  PRINT  "  ****  UERD/DINSRDC  (DNA).  «***" 

28  PRINT  "  ****  9  OCT,  1980  REVISION  TO  INCLUDE  ***«" 

30  PRINT  "  ****  EFFECTS  OF  BOTTOM  LAYER  ACCELERATION****" 

32  PRINT  "  *•**  BY  J.  CLARK  (CUA).  ***•" 

34  PRINT  "  ****  23  DEC  1987  REVISION  TO  INCLUDE  **«*" 

35  PRINT  "  ****  UPPER  BOUNDARY  DYNAMICS,  AFTERFLOW  «***" 

36  PRINT  "  *«**  AND  CLOSURE  MACH  NUMBER  **•*" 

37  PRINT  "  ****  BY  V.  CUSHING  (CAI).  *»**" 

38  PRINT  "  *HHm»m«**i»««mii»4»f*#iH»H»***»i«n 

39  RETURN 

40  WIDTH  LPRINT  100 

50  FF$  =  CHR$(  12) :  QU0TE$  =  CHR$(34):  CR$  =  CHR$(  13) 

60  DEF  FNNU$(I,J)  =  RIGHT$( SPACE$(J-1 )+STR$(I) , J) 

70  DEF  FNN$( I )  =  MID$( STR$( I ) ,2) 

80  "  KEY  7, "SAVE"+QUOTE$+"L.L"+QUOTE$+",A"+CR$ 

90  "  KEY  9,  "SYSTEM"+CR$ 

100  "  KEY  8.  "LPRINT  FF$"  +  CR$ 

110  DEFINT  I-N:  DEFSTR  S 

120  ;  DIM  R( 200) ,  Z(200,3),  T(200,3).  ZSD0T(200),  ZBD0T(200),  DELD0T(200), 
RMACH(200) 

122  N1  =  5:  N2  =  4:  SF0RM2  =  STRING$( N1 , "#" )  +  +  STRING$( N2 , "#" ) :  NN  = 

N1+N2+1 :  SF0RM1  =  "\"  +  SPACE$(NN-2)  +  "\" 

500  '  INITIALIZATION 

510  PATMOS  =  14.7*144:  RHO  =  2:  GEE  =  32.2:  CWAT  =  5000 

520  GAMMA  =  RHO*GEE :  AINT  =  2*RH0/PATM0S  *  < - >  PSI  IN  INTEGRAL  (TOTAL  IN 

LINE) 

530  RSTART  =  200 
540  AIMP  =  RHO*CWAT 

550  AK  r  144*23976:  AL  =  1.13:  AM  =  .0654*. 001:  AN  =  -.18  '  HI 

PARAMETERS 

554  GOSUB  20 

555  INPUT  "  OUTPUT  FILE  NAME  :  ",  SOUT 

556  INPUT  "  POUNDS  OF  HE  :  ",  W 

570  W3  =  W" .333 

574  INPUT  "  BURST  DEPTH  :  D 

580  PRINT  "  PROGRAM  IS  NOW  COMPUTING  MAXIMUM  CAVITY  RADIUS" 

590  NDZ  =  500 
1000  »  ENTRY  POINT 
1010  GOSUB  8000 

1020  PRINT  "  THE  MAXIMUM  CAVITY  RADIUS  IS  :  "  R(  0) 

1022  INPUT  "  #  RADIAL  DATA  POINTS  AS  N#:  OR  DR  AS  D#:  ",  S 

1024  IF  LEFT$( S, 1 )  =  "N"  THEN  DELRAD  =  R( 0)/VAL(MID$( S,2) ) 

1026  IF  LEFT$(S,1)  =  "D"  THEN  DELRAD  =  VAL(MID$( S,2) ) 

1030  PRINT  "  PROGRAM  IS  NOW  COMPUTING  VERTICAL  BOUNDS":  LIN 

CSRLIN :  GOSUB  9000 

1040  PRINT  "  PROGRAM  IS  NOW  COMPUTING  CLOSURE":  LIN  = 

LIN+1 :  GOSUB  10000 

1044  PRINT  "  PROGRAM  IS  NOW  WRITING  TO  FILE": 

GOSUB  15000 

1050  •  PLOT  CLOSURE  WAVE 
1054  PRINT  SPACE$( N1-2) ; 

1060  PRINT  USING  SF0RM1 :  "  R  " :  "Z-UPPER":  "Z-LCWER";  "Z-CLOSE"; 

"T-CLOSE":  "ZSDOT";  "ZBDOT" ;  "1/MACH" 

1070  FOR  I  =  1  TO  IMAX 

1080  PRINT  USING  SF0RM2:  R(I);  Z( 1, 1 ) ;  Z(I,2);  Z(I,3)j  T(I,3);  ZS  DOT(I); 

ZBDOT(I) ;  RMACH(I) 

1082  NEXT  I 
1100  SYSTEM 
1200  END 

2000  ' - - SUBROUTINE  LCWER(R,Z) 

2010  Z  =  0:  GOSUB  7000:  FMEM  =  F:  XMEM  =  Z 
2020  Z  =  2.5:  I FLAG  =  0:  ICOUNT  =  0 


AL  = 


.0654*. 001:  AN  = 


OUTPUT  FILE  NAME 
POUNDS  OF  HE 

BURST  DEPTH 


HBX-1 


L/UUU  1  ULii  114  •  •  i-' 

PROGRAM  IS  NOW  COMPUTING  MAXIMUM  CAVITY  RADIUS" 


THE  MAXIMUM  CAVITY  RADIUS  IS  :  "  R(  0) 

#  RADIAL  DATA  POINTS  AS  N#:  OR  DR  AS  Of:  ",  S 
"N"  THEN  DELRAD  =  R( 0)/VAL(MID$( S,2) ) 

"D"  THEN  DELRAD  =  VAL(MID$( S,2) ) 

PROGRAM  IS  NOW  COMPUTING  VERTICAL  BOUNDS":  LIN  = 

UB  9000 

PROGRAM  IS  NOW  COMPUTING  CLOSURE":  LIN  = 

B  10000 

PROGRAM  IS  NOW  WRITING  TO  FILE": 
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6050 

6060 

6090 

6100 

6500 

6510 

6520 

6530 

6540 

7000 

7010 

7020 

7030 

7040 

7050 

7060 

7070 

70&0 

6000 

8010 


GOSUB  7000:  X  =  Z 
GOSUB  4000:  Z  =  X 
IF  IFLAG  =  0  GOTO  2030 
RETURN 
* 

» -  SUBROUTINE  UPPER{R,Z) 

Z  =  0:  GOSUB  5000:  FMEM  =  F:  XMEM  =  Z  •  5000  REQUIRES  R,Z 
Z  =  .1:  IFLAG  =  0:  ICOUNT  =  0 
GOSUB  5000:  X  =  Z 
GOSUB  4000  :  Z  =  X 
IF  IFLAG  =  0  GOTO  3030 
GOSUB  6500:  T  =  TAR 
RETURN 
• 

« - SUBROUTINE  FOR  F >  0 

XNEW  =  X  -  F* ( X-XMEM ) / ( F -FMEM ) 

XDIF  =  ABS(XNEW-X) 

XMIN  =  .05 

IF  (XDIF  <  XMIN)  THEN  IFLAG  =  1 

XMEM  =  X 

X  =  XNEW 

FMEM  =  F 

RETURN 

t 

' - SUBROUTINE  UPP1(R,Z,F) 

GOSUB  6000 

F  =  -PREL  -  AK«((W3/RA)lAL) 

RETURN 

t 

•  - SUBROUTINE  PRESS( R,Z,PREL) 

GOSUB  6500 

THETAI  =  AM«W3#((W3/RI)lAN) 

PIO  =  AK#( ( W3/RI)~AL) 

PINC  =  PIO*EXP(-TT/THETAI) 

UI  =  ( PIO/AIMP)*(  THETAI/TIO  +  ( 1-THETAI/TIO) •(  EXP(-TI/THETAI)  )  ) 
'  INCLUDING  AFTERFLOW 
PHYD  =  GAMMA*Z 

PREL  =  -(PINC  +  PATMOS  +  PHYD):  UA  =  PREL/AIMP 

UZ  =  -UI*(D-Z)/RI  +  UA*( D+Z)/RA :  UR  =  UI*R/RI  +  UA*R/RA 

RETURN 

' _ _  T  SUBROUTINE 

RI  =  SClR(R*R  +  (D-Z)l2):  RA  =  SQR(  R*R  +  (D+Z)l2) 

TAR  =  RA/CWAT :  TIO  =  RI/CWAT:  TI  =  TAR  -  TIO 

RETURN 

» 

•  - SUBROUTINE  LOW  1  (  R, ZB, G3) 

GOSUB  6000 

FI  =  (D  +  Z)/RA 

FI2  =  (RA  -  2*D*F 1 )/RI 

F  =  -(  1  +  FI2*(AN*RA/RI  -  AN  - 1 ) ) »PINC/(CWAT«THETAI) 

F  =  F  -  AL*FI2*PINC/RI 

F  =  F  +  GAMMA *F1  -  AL*PREL/RA 

RETURN 

i 

•  -  SUBROUTINE  TANPT(R3) 

R  =  RSTART 

GOSUB  3000:  GOSUB  7000:  RMEM  =  R:  GMEM  =  F 
R  =  2*R 

GOSUB  3000:  GOSUB  7000:  IFLAG  =  0 

X  =  R:  XMEM  =  RMEM:  FMEM  =  GMEM:  GOSUB  4000:  R  =  X:  RMEM  =  XMEM: 

GMEM  =  FMEM 
IF  IFLAG  =  0  GOTO  8060 
R(  0)  =  R 
RETURN 

t 

•  - SUBROUTINE  B0UNDS(R(0)) 

R  =  0:  I  =  1:  FINAL  =  R(  0)  +  .5*DELRAD 

LOCATE  16.1:  PRINT  "RADIAL  COMPUTATION  #:  ":  GOSUB  20000 
'  PRINT  COUNTER 
GOSUB  3000 

Z(I,1)  =  Z:  T(I, 1 )  =  TAR 
GOSUB  2000 
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9060  Z ( 1,2)  =  Z:  T(I,2)  =  TAR:  R(I)  =  R 

9070  IF  Z  >  Z2MAX  THEM  Z2MAX  =  Z 

9090  R  =  R  +  DELRAD:  1=1+1 

9100  IF  R  <  FINAL  GOTO  9020 

9110  IMAX  =  1-1 

9120  RETURN 

9130  ' 

10000  » -  SUBROUTINE  CCURVE 

10010  FOR  I  =  1  TO  IMAX 
10020  DZ  =  Z2MAX/NDZ 
10030  R  =  R( I) 

10040  GOSUB  20000 
10060  ' 

10100  1 - SET  UP  BOTTOM  FACTORS 

10110  Z  =  Z( 1.2) 

10120  GOSUB  6000 

10130  TA  =  TAR:  TE  =  TI:  Til  =  TIO 
10140  THETA  =  THETAI 

10150  C0EF1  =  -PIO/AIMP :  C0NST1  =  THETA/TI1 :  C0EF2  =  1  -  C0NST1 

10160  C0EF3  =  C0EF1*THETA :  C0EF4  =  C0NST1  -  1:  C0NST2  =  EXP(-TE/THETA) 

10170  FACTVI  =  (D-Z)/RI:  FACTVA  =  (D+Z)/RA 

10180  FACTRI  =  R/RI :  FACTRA  =  R/RA 

10190  UAO  =  UA 

10210  Z  =  0:  GOSUB  6000:  UZMEM  =  UZ 

10220  TOTAL  =  0:  ZMEM  =  0:  ZB  =  0:  ZS  =  0:  ZSMEM  =  0:  TTMEM  =  TAR: 
IFLAG2  =  0 

10230  GOSUB  10500:  FMEM  =  F 
10240  Z  =  DZ:  ICOUNT  =  0 
10250  ' 

10300  1  -  CLOSURE  COMPUTATION 

10310  GOSUB  11000 

10320  IF  I FLAG 2  =  1  GOTO  10360 

10330  IF  ZB  <  ZS  THEN  IFLAG  =  0:  IFLAG2  =  1:  GOSUB  10500:  FMEM  =  F: 
XMEM  =  Z:  Z  =  Z  -  DZ :  GOTO  10310  '  CHANGE  TO  F  ROUTINE 
FOR  CONVERGENCE 
10340  Z  =  Z  +  DZ 
10350  GOTO  10300 
10360  GOSUB  10500 
10370  X  =  Z:  GOSUB  4000:  Z  =  X 
IO380  IF  IFLAG  =  0  THEN  ZMEM  =  XMEM:  GOTO  10300 
10390  GOSUB  11000 

10400  IF  I  =  1  THEN  RMACH(I)  =  0:  GOTO  10420 

10410  DS  =  SQR(  ( R-RLAST )_2  +  (Z-ZLASTT2  ):  DT  =  TT  -  TLAST :  RMACH(I) 
DT*CWAT/DS 

10420  RLAST  =  R:  TLAST  =  TT:  ZLAST  =  Z 

10430  Z( I, 3)  =  Z:  T(I,3)  =  TT:  ZSDOT(I)  =  ZSDOT:  ZBDOT(I)  =  ZBDOT : 
DELDOT(I)  =  ZSDOT  -  ZBDOT 


10480  NEXT  I 
10490  RETURN 

10500  ' - SET  UP  F  FOR  CLOSURE 

10510  F  =  ZS  -  ZB 
10520  RETURN 
10530  ' 

11000  ' - START  LOOP  WITH  SURFACE  MOTION 

11010  GOSUB  6000 


11020  IF  IFLAG2  =  1  THEN  DZ  =  Z  -  ZMEM 

11030  TOTAL  =  TOTAL  +  .5*(UZ  +  UZMEM )*DZ:  UZMEM  =  UZ 

11040  TF  =  AINT*( Z*UZ  -  TOTAL)  '  TIME  OF  FLIGHT  TF  (SEC)  AS  FUNCTION  OF 
SPALL  THICKNESS  Z  (FEET) 

11050  ZR  =  UR*TF 

11060  ZS  =  TF*UZ  +  .5*GEE»TF»TF  '  SURFACE  DISPLACEMENT  —  EQS.  14  AND 
18  REFERENCE  15 

11070  TT  =  TAR  +  TF  '  TT  =  TOTAL  OR  ACTUAL  TIME 

IIO80  DELT  =  TT-TTMEM 

11084  IF  ABS(DELT)  <  .000003  GOTO  11100 

11090  ZSDOT  =  ( ZS-ZSMEM) /DELT 

11100  ZSMEM  =  ZS:  ZMEM  =  Z:  TTMEM  =  TT 

11110  IF  TT  <  TA  THEN  Z  =  Z  +  DZ :  GOTO  11000 

11200  '  -  BEGIN  BOTTOM  MOTION  CALCULATION 

11210  T  =  TT-TI1 

11220  VI  =  C0EF1*(  C0EF2*EXP( -T/THETA)  +  C0NST1  ) 

11230  ZI  =  C0EF3*(  (  EXP( -T/THETA)  -  C0NST2  )*C0EF4  +  (T-TE)/TI1  ) 
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'  REFERENCE  11,  EQ.  23  —  NEGATIVE  (UPWARD)  DISPLACEMENT 
11240  T  =  TT-TA 
112S0  VA  =  UA0*( 1  +  T/TA) 

11260  ZA  =  UA0*T*( 1  +  .5*T/TA)  '  EQ.  24 

11270  ZBR  =  ZI*FACTRI  +  ZA*FACTRA:  ZBRDOT  =  VI*FACTRI  +  VA*FACTRA 
11260  ZB  =  ZI*FACTVI  +  ZA*FACTVA :  ZBDOT  =  VI»FACTVI  +  VA*FACTVA 
11300  RETURN 
11310  ' 

15000  ' - SUBROUTINE  DFILE 

15010  OPEN  "0",  #1,  SOUT 

15011  PRINT  #1,  SPACE$(N1-2); 

15012  PRINT  #1,  USING  SF0RM1;  "  R  nZ-UPPERn;  "Z-LCWER":  "Z-CLOSE";  "T- 

CLOSE":  "ZSDOT";  "ZBDOT":  "1/MACH" 

15020  FOR  I  =  1  TO  IMAX 

15030  PRINT  #1,  USING  SF0RM2;  R(I):  Z(I,1);  Z(I,2);  Z(I,3);  T(I, 3) ;  (I); 

ZSDOT(I) ;  ZBDOT(I) ;  RMACH(I) 

15044  NEXT  I 
15046  CLOSE  # 1 
15050  RETURN 
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